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Energy extremum principle for charged black holes 

Scott FraseiQ and Shaker Von Price FunkhouseiHI 
Department of Physics, California Polytechnic State University, San Luis Obispo, California 93407, USA 

For a set of N asymptotically flat black holes with arbitrary charges and masses, all initially at rest 
and well separated, we prove the following extremum principle: the extremal charge configuration 
{\qi\ = nii for each black hole) can be derived by extremizing the total energy, for variations of the 
black hole apparent horizon areas, at fixed charges and fixed Euclidean separations. We prove this 
result through second order in an expansion in the inverse separations. If all charges have the same 
sign, this result is a variational principle that reinterprets the static equilibrium of the Majumdar- 
Papapetrou-Hartle-Hawking solution as an extremum of total energy, rather than as a balance of 
forces; this result augments a list of related variational principles for other static black holes, and is 
consistent with the independently known Bogomornyi-Prasad-Sommerfield (BPS) energy minimum. 

PACS numbers: 04.70.Bw, 04.40.Nr 


I. INTRODUCTION 

A system of multiple black holes is not static in gen¬ 
eral, but the system can be kept static by a strut (a con¬ 
ical singularity) on the axis between each pair of black 
holes. Each strut is interpreted as providing a force that 
prevents the black holes from moving (see [1] and refer¬ 
ences therein). Remarkably, without any struts present, 
there is also a well-known solution [5] that describes a 
set of N charged black holes in static equilibrium, re¬ 
gardless of their mutual separations. A key feature of 
this static solution is that the black holes’ individual 
charges qi all have the same sign, and are related to 
their individual masses m,- by the condition (in units with 
G = c = dTTCo = 1) 

l^il = m*, l<i<N. (1) 

The asymptotically flat static geometry satisfying Q 
was found by Hartle and Hawking [2], who extended the 
Majumdar-Papapetrou solutions [all!. The interpreta¬ 
tion and significance of the condition 0 is of central 
interest in this paper. The traditional interpretation of 
the Hartle-Hawking solution [2] appeals to forces: the 
static equilibrium is attributed to the exact balance be¬ 
tween the black holes’ pairwise gravitational attraction 
and electric repulsion. However, the condition Q refers 
to each black hole, not to a balance of forces between 
black holes. Thus, 0 is a stronger condition than a New¬ 
tonian balance of forces, which for two particles requires 
qiq 2 = TO 1 TO 2 . This suggests that ([^ can be interpreted 
without using forces. 

Such an alternative interpretation is the purpose of 
this paper. We show that the condition ([^ can be de¬ 
rived from an extremum of the system’s total energy E. 
This fits naturally into the framework of general relativ¬ 
ity, where gravity and total energy are determined by 
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spacetime geometry, and gravity is fundamentally not 
treated as a force. 

The condition 0 describes extremal charge, since a 
nonrotating black hole satisfies a charge-mass inequality, 
|g| < TO. Thus, our results can be briefly summarized as: 
a set of extremally charged black holes extremizes the 
total energy. More precisely, we consider a conformally 
flat spatial geometry in the form introduced by Brill and 
Lindquist [5| , and we prove the following extremum prin¬ 
ciple for well-separated black holes: 

For N charged black holes, initially at rest, each 
with \ qi\ < rrii and apparent horizon area Ai, the 
extremal charge condition |9i| = to^ follows from 
extremizing the total energy: dE/dAi = 0 at 
fixed charges qi and Euclidean separations Vij. (2) 

We prove the extremum principle (1^ as an expansion in 
the inverse separation distances (l/r-y), through second 
order, which is where relativistic post-Newtonian contri¬ 
butions first appear. The particular quantities that are 
varied or held fixed in this energy extremum are moti¬ 
vated by the first law of black hole mechanics, as we will 
illustrate when we prove ([^. 

If all of the black hole charges in 0 have the same 
sign, then ([^ is a variational principle for static black 
holes: it identifies a static black hole configuration as an 
extremum of total energy, within a family of black holes 
that are initially at rest. That is, a configuration with 
extremal energy remains at rest, while all other config¬ 
urations evolve dynamically. Such variational principles 
have been proved for a single uncharged black hole |6] , a 
single black hole in Einstein-Yang-Mills theory [7|, and 
pairs of mirror-symmetric black holes in the Randall- 
Sundrum braneworld models |S]. Our extremum prin¬ 
ciple ([^ adds the case of charged black holes to these 
earlier variational principles ISHH] and identifies the static 
Hartle-Hawking solution as its energy extremum. In con¬ 
trast to (in, in this paper, the number N of black holes 
that we consider is arbitrary, and our energy extremum 
varies (rather than holds fixed) the black hole areas. 
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We also show that the energy extremum in ([^ is an en¬ 
ergy minimum. This result is consistent with the earlier 
result of Gibbons and Hull [5] , but our methods are differ¬ 
ent. In [^, it was found that the Majumdar-Papapetrou 
solution saturates the lower bound of a supersymmet¬ 
ric Bogomornyi-Prasad-Sommerfield (BPS) inequality, 
E > IQI, for global quantities: the total energy E and 
charge Q (with zero magnetic charge in the context of 
this paper). In contrast to [5], in this paper, we min¬ 
imize the energy E(Ai,qi,rij) as a function of several 
variables, and we do not use supersymmetry; we also ob¬ 
tain additional extremum conditions (m^ = |gi|) on the 
individual masses and charges. 

It is worth noting that, although black holes with ex¬ 
tremal charge are not often considered in astrophysi- 
cal applications, they have attracted significant theoret¬ 
ical interest. For example, they represent stable ground 
states for black holes in supersymmetric theories [10] and 
they have been used to investigate the nature of entropy 
in black hole thermodynamics m- 

This paper is organized as follows. In Sec. [TTj we re¬ 
view the necessary geometry. In Sec. |III[ we prove the 
extremum principle © for well-separated black holes, 
through second order, and we specify the physical con¬ 
ditions under which the separations Vij are sufficiently 
large. In Sec. EYl we show that the energy extremum 
is a minimum, and verify that it agrees with the BPS 
bound. We conclude in Sec. |Vj Throughout this paper, 
we work in four spacetime dimensions and use geometric 
units with G = c = 47reo = 1- 

Much of our analysis does not require that the black 
hole charges have the same sign. We only need to refer to 
same-sign charges in the following contexts: in a higher- 
order analysis (Sec. |IIIE ), in comparisons to the BPS 
bound (Sec. |IV| ), and in applications of the extremum 
principle ([^ as a variational principle that reproduces 
the static Hartle-Hawking solution (as described above). 


II. GEOMETRY 

In this section, we review the geometry mis] that we 
use in this paper. A system of N charged black holes, all 
initially at rest, is described by its instantaneous spatial 
geometry. The appropriate area Ai of a black hole is 
that of its apparent horizon, which is determined by the 
spatial geometry alone (unlike the event horizon, which 
is a global spacetime property). The apparent horizon 
generally lies inside the event horizon, and coincides with 
it for a static or stationary spacetime. 

We use a conformally flat geometry, in the form in¬ 
troduced by Brill and Lindquist [5]. This geometry con¬ 
tains no conical singularities (struts) to prevent the black 
holes from moving. The extremal black holes consid¬ 
ered by Hartle and Hawking [2] remain eternally static. 
The nonextremal black holes considered by Brill and 
Lindquist [S] are initially at rest, and evolve dynamically 
thereafter [H For both cases, the instantaneous 



FIG. 1. Illustration of three nonextremal black holes in the 
spherical approximation. Black holes 1 and 2 are separated 
by distance ri 2 . Black hole i has position iy and apparent 
horizon radius Ri. Dots denote the punctures described in 
the text. 


spatial geometry, exterior to all N black holes, is 

ds^ = [dx^ + dy^ + dz^) , (3a) 


/ 




(3b) 


Figure illustrates the setup. The Euclidean vector r 
locates any point outside the black holes. The vector 
r) locates black hole i. We let denote a Euclidean 
separation distance between black holes i and j. 


Tij = \ri-rj\. 


(4) 


The parameters ai and fdi are non-negative constants, 
and are related to physical quantities (mass, charge, en¬ 
ergy) as follows [21 [5] . Black hole i has individual mass 
(rest energy) nii and charge qi, 




q^= j3i-ai + y^^ 




(5) 

( 6 ) 


Each pair {mi,qi) satisfies the black hole charge-mass 
inequality, \qi\ < mi. The system of N black holes has 
total energy E and interaction energy Ejnt, 

N 

E = ^(ai + /3i), (7) 

= = ( 8 ) 

i—1 2—1 j^i 

The expressions for rui and qi^ in ^ and ®> , can be 
solved for the parameters ai and /?i. For large Vij, the 
results to first order in l/r^ are [5] 


di — 


(mi - qi) 


{rrii + qi) 


^ _ 1 [mj -h qj) 
2 ^ fij 

^ _ 1 [mj - qj) 


2 Tij 


(9a) 

(9b) 
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Using , for large , the interaction energy Q can be 
written to first order in 1 /r^ as [5] 


N 


i=l j>i 


{m^rrij - q.,qj) 


( 10 ) 


This has the expected form: it is the sum of the pairwise 
gravitational and electrostatic potential energies. 

For the static Hartle-Hawking solution [5] , only one set 
of the parameters Ui and /3i are nonzero, 

(ai=0,/3i>0) or {/3^ = 0,a^ > 0). (11) 

In terms of mass and charge, the two cases in 0 are 
equivalent to, respectively [2] 


qi = rrii = Pi > 0 


or 


- qi = mi= ai> Q. ( 12 ) 


Thus, the static condition, (111 and (12), can be summa¬ 


rized as an extremal condition for same-sign charges, 

\qi\ = iTii, with all qi of like sign. (13) 

The black hole charges are nonextremal, if the pa¬ 
rameters at and Pi do not satisfy the static condition 
(11|-(13). Near a nonextremal black hole i, we may 
transform from Cartesian coordinates {x,y,z) in ^ to 
spherical coordinates centered at f(. As found 

in [^, for sufficiently large separations r^, the appar¬ 
ent horizon of black hole i can be treated as a sphere, 
r{9, (j>) = Ri = constant. We will refer to this as the 
spherical approximation, which is illustrated in Fig. 
The radius Ri is given by [5] 






(14) 




In the spherical approximation, the metric function / in 
^ evaluated on the surface r = Ri is [5] 


/H'+i+gilh+t+gli- (-) 


As we see in Sec. Ill the results (14) and (15) will only 


hold to up to an appropriate order, when expanded in 
powers of the inverse separations (l/r^). 

We end this subsection by reviewing some technical 
features [H [5] of the geometry and topology, for both 
nonextremal black holes and extremal black holes. The 
spatial geometry has the topology of with N points 
(‘punctures’) removed. Each puncture is located by the 
vector -Ti. In the flat background metric, the distance 
between two punctures is the Euclidean distance r^-. 

In the full conformal geometry, for a nonextremal black 
hole, each puncture represents the spatial infinity of an 
asymptotically flat region (‘sheet’), which is hidden be¬ 
hind the black hole’s apparent horizon, as viewed in the 


common asymptotic region (sheet) exterior to all N black 
holes, where the total energy E is dehned. The topology 
referred to above is equivalent to the topology of -I- I 
connected sheets. Each black hole’s individual mass rrii 
is the Arnowitt-Deser-Misner (ADM) mass of its hidden 
asymptotic region (its individual sheet). 

In contrast, for an extremal black hole, a careful co¬ 
ordinate analysis P] reveals that the puncture at fy rep¬ 
resents the apparent horizon, with finite area. In this 
case, the parameter values 0 cannot be used to eval¬ 
uate the apparent horizon radius (14). As shown in P], 
for extremal charge, |r — rij —?► 0 is a simple coordinate 
singularity, and the area of the apparent horizon in this 
limit is nonzero, Ai = AKurii^ = Ai^qp. This coordinate 
issue will pose no complications in our proof of the ex¬ 


tremum principle in Sec. Ill B where we will perform all 


calculations in the nonextremal regime, and then express 
the energy E as function of the finite areas Ai. These 
areas reproduce the correct values Ai = AttitiP = ^irqp 
when we take the extremal limit in terms of charges and 
masses (|gi| —>■ rrii), which is a coordinate-independent 
limit. 


III. ENERGY EXTREMUM PRINCIPLE 
A. Physical motivation 

In this subsection, we illustrate how the first law of 
black hole mechanics motivates our extremum principle 
([^. This also demonstrates the basic extremum proce¬ 
dure that we use in our proof of the extremum principle 


(Secs. IllB THE). 


For clarity, we begin with the case of a single black 
hole, and then generalize it to the multi-black hole 
case. A single static charged black hole (the Reissner- 
Nordstrom solution) has energy E = m and charge 
|g| < m. Variations conserve energy via the first law 
of black hole mechanics [14] , 


SE = —5A + $(5(7. 
Stt 


(16) 


Here A, k, and $ are the black hole’s surface area, surface 
gravity, and electric potential, respectively. In the case 
of extremal charge {\q\ = m), it is well known that At = 0, 
so the first law reduces to 


6E = $(5q. 


(17) 


The key point of 0 is the following: if |( 7 | = m, then 
variations that hold q constant {6q = 0) extremize the 
energy {SE = 0), and this means dE/dA = 0, since E is 
a function of A and q. Conversely, if q is held constant 
in ( |16| ), then an energy extremum requires dE/dA = 0, 
and this reduces to |( 7 | = m, as follows. The total energy 
of the Reissner-Nordstrom solution is 


( 18 ) 
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Evaluating dEjdA = 0 for (18) yields |g| = yjAjAir, for 
which E = m = \/AjAir = |q| and hence |g| = m. In 
summary, the extremum dEjdA = 0 occurs if and only 
if the extremal charge property holds, |g| = m. 

As already noted, the vanishing coefficient of 5A in the 
first law ( |17[ ) motivates why, for an arbitrarily charged 
black hole, we extremize the energy (18) by varying the 
area A and holding the charge q constant. This general¬ 
izes to the multi-black hole case, as follows. The expected 
first law for the static Hartle-Hawking solution [5] is 


C. Proof through first order 

We begin by calculating the apparent horizon area Ai 
of each black hole. As in the original work of Brill and 
Lindquist [S], through first order in l/r^j, we use the 
spherical approximation (see Sec. for which the ap¬ 
parent horizon is the surface r = Ri = constant in spher¬ 
ical coordinates {r,9,(j)) centered at ry. Then the differ¬ 
ential area is 


dAi = f^Ri^ sin^ 6 dO dcj, (20) 


N 

= (19) 

i=l 

where E is the total energy. For arbitrarily charged black 
holes, in analogy with the single-black hole case, we ex¬ 
tremize E by varying the black hole areas Ai while hold¬ 
ing constant all other quantities: the charges qi (which 
appear in the first law) and also the black hole separa¬ 
tions rij (which are internal parameters of the geometry). 
This motivates our extremum principle ©• 

Note that the single black hole considered above re¬ 
mains at rest, so it illustrates our energy extremum pro¬ 
cedure, but not any subsequent dynamics. In this paper, 
the initially static multi-black hole configurations that 
are nonextremal will not remain static; their dynamics 
were studied in [12]. 


B. Outline of proof 


where Ri and / are given by (14) and (15), respectively. 
Integrating then yields the area formula: 

A, = AT:R,^f. ( 21 ) 


We now evaluate this using (14) and (15). Expanding 
through first order in l/r^ gives 





( 22 ) 


with rrii given by ([^ as a function of the other parameters 
{ai, /3i,rij). We can express the area Ai completely in 
terms of the mass rui and charge qi, by substituting ai 
and in ([^. Then simplifying yields 



= mi + \lmj^ - qj^ 


X 





(23) 


In Secs, jlll CHIII E| below, we will prove the extremum 
principle (plas an expansion in the inverse separations, 
l/cy. Zeroth order (l/r^ —)■ 0) corresponds to —>■ oo, 
and treats each black hole in isolation from the other 
black holes. Through first order, we retain all terms lin¬ 
ear in 1/rij, which characterize effectively Newtonian in¬ 
teractions among the black holes. Through second order, 
we retain all terms quadratic in 1/, which characterize 
post-Newtonian relativistic interactions. 

All of the terms in our expansions will be dimensionless 
ratios of black hole properties {ai, j3i, rrii, qi) divided by 
the separations rij. These ratios are small for sufficiently 
large separations. In particular, we will show that the 
extremum principle holds for sufficiently well separated 
black holes, as specified by conditions on rij that take the 
form Ui > 0, where Ui can be interpreted as the effective 
gravitational potential experienced by black hole i. 

At each order, the main steps are to calculate each 
black hole’s apparent horizon area Ai, and then to ob¬ 
tain the energy E{Ai,qi,rij) as a function of the areas, 
charges, and separations. We then extremize the en¬ 
ergy as dEjdAi = 0, and from this obtain the condition 
\qi\ = rrii, which proves the extremum principle. 


When we expand this expression, the terms that are first 
order in l/r^- exactly cancel. At first order, we neglect 
the terms quadratic in l/cj^, and obtain 


^ = rrii + - qj^. 


(24) 


Solving this for the mass rrii yields 
TT f Ai 



(25) 


Due to the cancellation that occurred to obtain (24), the 


mass (25) contains neither first-order terms (proportional 
to 1 jrij) nor any parameters related to the other black 
holes {j ^ i). This first-order result is therefore the same 
as the zeroth-order result, and it has the same form as 


the energy (18) for a single black hole. 


We now evaluate the total energy E, which by (^8| and 


(10) is, to first order in l/r^, 


N N 

i—1 j>i 


{rriimj - qiq^) 


( 26 ) 


2=1 





















5 


The zeroth-order term in (261 is the sum of the rest en¬ 


ergies, and the first-order terms are the pairwise gravi¬ 
tational and electrostatic potential energies; this will be 
useful in interpreting the quantity Ui below. 


We now take each mass rm in (26) to be given by (251, 


which expresses the energy (26) in the functional form 


E{Ai,qi,rij). Then holding constant the charges qi and 
separations rij, we find 


where 


dE _ dm^ 
dA ~ dA,, 




(27) 


(28) 




We now extremize the total energy E by requiring 
dE/dAi — 0, which by (27) is equivalent to 




(29) 


For sufficiently large separations, Ui ^ 0, as seen by (28). 
In this case, the extremum (29) requires 


drrii 


= 0 . 


(30) 


Evaluating this for (p^ reduces to 


For this charge valui 5) is 



(31) 


(32) 


Thus, since the extremum ( |29[ ) has produced the N de¬ 
sired conditions, \qi\ = rui, we have proved the extremum 
principle ([^ through first order. 

As noted above, our proof requires Ui ^ 0. As seen 
by (p^, this condition is always satisfied at zeroth or¬ 
der (T/rij — )■ 0 and Ui = 1), and it is satisfied at first 
order for sufficiently large separations compared to 
the masses rrii. In particular, the continuity of our per¬ 
turbative approach with the zeroth-order limit {Ui = 1) 
requires 0 < Ui < 1. 

We end this subsection by interpreting the quantity Ui. 
Note that Ui can be obtained from the energy (26) in the 


form E{mi,qi,rij), along with (27) and the chain rule. 


U, = 


dE 

dnii 


(33) 


This expression motivates the interpretation of Ui as the 
effective gravitational potential experienced by black hole 
i, including its rest energy. This interpretation is justi¬ 
fied by the derivation of (27)-(28) from pM, which makes 


clear that: (i) the zeroth order term in ( |28| ) is the rest en¬ 
ergy, normalized per unit mass; (ii) the first-order sum in 
(28) is the Newtonian gravitational potential; (iii) there 
is no electric potential contribution to Ui since (27|-(28) 


are obtained by holding the charges constant in (26). 


z 



FIG. 2. Illustration of two nonextremal black holes in the 
second-order analysis, showing useful distances and coordi¬ 
nates {r,9). The dots have Euclidean separation rij = \U—ri\. 
All black holes (including others, not shown) are collinear on 
the 2 axis, and have azimuthal symmetry about the a axis. 


D. Proof through second order 


In this section, we prove the extremum principle (pj) 
through second order. Our first task will be to evaluate 
each apparent horizon area Ai. This is a longer calcu¬ 
lation than in our first-order proof (Sec. IIIC), since at 
second order, we must depart from the spherical approx¬ 
imation. Figure p] illustrates the setup. 

In general, each black hole’s apparent horizon is non- 
spherical, due to its interactions with the other black 
holes. Thus, in spherical coordinates {r,0,(j)) centered 
at fi, the apparent horizon of black hole i is generally 
a surface r{9,(j)). To simplify the analysis, we will take 
all black holes to be collinear (aligned on the z axis). 
The resulting axisymmetry simplifies each nonspherical 
horizon r{9,(j)) to the azimuthally symmetric form r{9). 
Thus, as illustrated in Fig. pj we take all N black hole 
positions fi to lie on the z axis. Near black hole i, we in¬ 
troduce spherical coordinates {r,9,(j)) centered at f, and 
consider another black hole labeled by j ^ i. This gives 


\r-r,\ =r, 


\f-fj \ = Jr'^ - cos 9 + Vij'^. 


(34a) 

(34b) 


To calculate the area Ai, we only need to consider points 
{r,9) on the apparent horizon, so r < r^. In this case. 


polynomials Pk, 


the distance (34b) has a standard expansion in Legendre 


r - Ti 




T h \ T< 

k>0 


( 35 ) 
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The form ( |35[ ) is very suitable for our purposes, since it 
is an expansion in the inverse separations (1/rij). Our 
first-order proof (Sec. Ill C I, like the original work of Brill 
and Lindquist [5], used the summation in ( [3^ through 
k = 0 and Pq = 1- Here, our second-order proof proceeds 
through k = 1 and Pi = cos 6*. 

^1) and ([^ to evaluate the metric func- 
We retain all terms 


We now use 
tion / in ([^ near blac 


:k hole i. 

through second order in Ijrij. This gives 


f{r,0) = l + 


Oii -j- [3i OCij3i 


E 


1 r cos 9 


E Oij Pk 

T* ' ■ ^ T* ' 7 

. , . ' LI j , . ' ik 


'"ij fij 


X [aj + Pj + 


^iPj T PiOlj 


• (36) 


This contains all of the second-order terms from expand¬ 


ing the spherical approximation (15), and generalizes it 
by including nonspherical terms, proportional to cos 9. 

Each black hole’s apparent horizon r{9) is an extremal 
surface, so it is determined by extremizing the area [S], 

pTT 

Ai = 2TT / d9 sm9r\/f^, (37) 

Jo 


where r = dr/d9. Extremizing (37) gives 


r + r 
1 -I- 


= —r ( cot 9 -\- 


2def 

f 


2rdrf 

f 


(38) 

As in jS], we will solve (38) by expanding in Legendre 
polynomials P„ as basis functions. This function-basis 
expansion for r{9) is distinct from our previous multipole 
expansion (35). Thus, we take 


^(^) = E CniPnjcosd), 


(39) 


n>0 


where Cni are constant coefficients and i labels the black 


hole. Our first-order proof (Sec. IIIC), like the original 


through n = 1, 


r{9) = Coi + CiiCos9, 


work of Brill and Lindquist [5] , used the expansion (39) 
through n = 0, which is the spherical approximation (14). 
At second order, we use 


(40) 


with Pg = 1 and Pi = cos 6*. To determine the coeffi¬ 


cients Cni, we evaluate the differential equation (38) for 


r{9) using (40) and the metric function / in (36). After 
expanding all quantities to second order in ijrij, (38) 
then takes the form 


FqPq + PiPi(cos0) — 0, 


(41) 


where Pg and Pi are functions (which for brevity we do 
not list here) that are independent of 9, but involve the 
other parameters (a^. Pi, rij). Since (41) is an expansion 


in the basis functions Pn, both Pg and Pi must vanish. 
Solving the equation Pg = 0 reduces to 


Coi Ri, 


(42) 


where Ri is the radius in the spherical approximation 
(14), which is here to be expanded through second order 


in 1/rij. Solving the equation Pi = 0 yields 


Cu = - 


^iPi 


4 -\/OiiPi + ai - 

2i/ aiPi{aj 


- Pj) + aiPj 


PiCXj 


(43) 


Thus, the coefficients Cu are second-order quantities. In 
[5], the coefficients Cni were determined numerically for 
uncharged black holes, and \Cni\ was found to decrease as 
n increases. Our results, (42) and ( |4^ , show analytically 
that here \Cni\ decreases as n increases due to its leading 
power in the inverse separations, (l/r^) 


2n 


For the coefficients Cni in (42) and (^43^ , the horizon 


surface (40) approximates a prolate ellipsoid with one 


focus at r = 0. This can be seen by linearizing the surface 
function re(0) of a prolate ellipsoid in the eccentricity 
(e< 1), 


m = 


n(l-e^) 

1+6 COS 6 


a — ae cos t 


(44) 


The linearized form matches the horizon (40), so the 


horizon approximates an ellipsoid with semimajor axis 
a = C^i = Ri and eccentricity e = Cu/Ri, as illustrated 
in Fig. [§ 


We now evaluate the horizon area Ai in (37) for the 


surface r{9) in ([40|, using the metric function / in (36). 


Expanding through second order in l/r^j, we find 






■\J Oli^i 


E 


4 \ ^ rij 


(45) 


with rrii given by ([^ as a function of the other parameters 
{ai, Pi,rij). The result (45) is the first-order area (1^, 


plus second-order corrections. 

Our next step is to obtain the energy E{Ai,qi,rij) as 
a function of the areas, charges, and separations. We 
evaluate the energy P by a different method than in our 
first-order proof (Sec. IIIC), where we used the Newto¬ 


nian form (26) of the energy, and expressed the masses 


rrii in terms of the areas Ai. At second order, our setup 
is post-Newtonian, and we find it convenient to express 
the parameters and Pi in terms of the areas Ai and 
charges qi, and then evaluate the energy ([^ in the form 

F = J2ii^i + A)- 

Thus, we must solve the area equation ( |4^ and the 
charge definition (§ for the parameters (ai, pi) as func¬ 
tions of the areas and charges {Ai,qi). For convenience 
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in summarizing our results below, we define the follow¬ 
ing functions Qi) of the areas and charges. We 

define the quantity as 




(46) 


which is our mass result (251 through Hrst order. We 

(47) 


define the two dimensionless quantities (Mi, Qi) as 
^^3 n 


11 11 


We now solve for the parameters (ai, Pi) using a pertur¬ 
bation approach: we write ai and Pi as their leading- 
order values plus second-order corrections (Pi, Pi). 
We rephrase the leading-order values (|^ using the defi¬ 
nitions (|46|) and (|47|. This gives 


(^J'i - li) 
2 

(Mz + li) 


1 - 


1 - 


(M^ + Q.) 

2 

(M, - Q.) 


+ Oii, 

+ Pi- 


(48a) 

(48b) 


It remains to solve for Pi and Pi as functions of the areas 
and charges (Ai,qi). To this end, we insert (1^) into the 
area result ( [4^ and into the charge dehnition (^. After 
expanding through second order in we then find 

that (45) and reduce to, respectively, 


K+ai +K^Pi = Ei, Pi-ai = Si. (49) 
We will summarize the coefficients K± and source terms 


(ei, (5i) below. Solving the linear equations (491 yields the 
second-order corrections. 


£i - K_5,, 


ft = 


Ei -|- K^6i 


K++K_' " K+ + K_ 

The coefficients K± and source terms (Ei,Si) are 

_ j -^i o q- 

— \ ft Qi I 

V 47r 

— Qi^i 


(50) 


(51) 

(52) 

(53) 


For convenience, we have let Si and 7) denote the follow¬ 
ing dimensionless second-order quantities. 


ft = 


Mp' Qp (l^jfJ-k — ijlk) 

4 Ararik 

k^j ^ ^ 

(fj-jQk - qJ^kk) 


7 ; = ft^-EE 

j^i kiij 


AVijTjk 


(54a) 

(54b) 


Two simple results follow from using the second-order 
solutions for ai and Pi. First, when we evaluate the mass 


rrii in ([^ using (48)-(54), a significant number of terms 
cancel, and we find 


rrii = 


(55) 


with Hi defined in (46). Interestingly, this is the same as 
the mass function (25) through first order. Henceforth, 


we wri te m , in place of A second result is that, after 
using (50)-(54), the sum Pi -|- ft reduces to 


«! + ft = TOift - gift. 


(56) 


We now evaluate the total energy E = + ft)) 

given by ([^. Using (48) and (56) gives 


N r 


^ = E 


rrii - ^(rritMi - qiQi) + nnSi - qiTi 


. (57) 


The zeroth-order and first-order terms in (57) reproduce 


our previous result (26), which is the sum of the rest 


energies and the pairwise Newtonian and electrostatic 
potential energies. The second-order terms in (57), pro¬ 


portional to ft and ft, are relativistic post-Newtonian 
corrections. 

The energy ( [A7| ) can be regarded as a function 
E(Ai,qi,rij) of the areas, charges, and separations. This 
follows from using the definitions of the quantities shown 
(rrii. Mi, Qi,Si,7i). Holding constant the charges qi and 
separations rij, we then find 


dE 


= ft 


duii 


(58) 


where 


ft = 1 - A4, + - (A4,2 + 

l 3 (rrijiTik — qjqk) 

njTjk 

(rrijmk - qjqk) 


Mi Mr 

ly y 

y. TikTjk 

3 k^t,j 


(59) 


In (59), the leading term (1 —AI^) is our first-order result 
(28) for ft. The additional terms in (59), which are sec¬ 


ond order, are relativistic post-Newtonian contributions. 

The remaining steps in our proof are now essentially 
the same as (29)-(32) in our hrst-order proof (Sec. HI C). 
We extremize the total energy E in (|57[) by requiring 


dEjdAi = 0, which by (58) is equivalent to 


ft^-0 


(60) 


For sufficiently large separations, ft ft 0. In this case, 
the extremum (60) requires dtrii/dAi = 0. Evaluating 


this for (46) reduces to \qi\ = Ai/Arr. For this charge 
value, the mass (46) is rui = Ai/An = \qi\. 

Thus, since the extremum (60) has produced the N 


desired conditions, {qp = rrii, we have proved the ex¬ 
tremum principle ([^ through second order. As noted 




























































above, our proof requires Ui ^ 0. This condition is al¬ 
ways satisfied at zeroth order (for which Ui = 1), and 
it is satisfied through second order for sufficiently large 
separations Vij compared to the masses and charges, as 
seen by (59), and by (62) below. As in the first-order 


case, continuity of our perturbative approach with the 
zeroth-order limit {Ui = 1) requires 0 < 17^ < 1. 

We end this subsection by revisiting the interpretation 
of Ui- As in the hrst-order case (Sec. IIIC), note that 
Ui can be obtained from the energy w in the form 
E{mi,qi,rij), along with (58) and the chain rule, 


U,= 


dE 

dnii 


(61) 


As in the first-order case, this motivates the interpreta¬ 
tion of Ui as the effective gravitational potential experi¬ 
enced by black hole i, including its rest energy. In this in¬ 
terpretation, it is not surprising that the gravitational po¬ 
tential (59) contains quadratic charge-dependent terms; 


this is because nonzero electric stress-energy (a quadratic 
function of the electric field) contributes to curving the 
geometry, and so contributes to the gravitational held. 
As an explicit example, for two black holes {N = 2), 


C/i = 1 - — 
ri2 


m 2 TO2^-I-92^ + 2(toiTO2 — 9192) 


2ri22 


, (62) 


and U 2 is similarly given by interchanging all subscripts 
(1 •(->■ 2). If we interpret Ui as an effective gravitational 
potential, then the terms proportional to toiTO 2 — qiq 2 
(which refer to black hole 1 itself) represent the nonlinear 
gravitational coupling in general relativity between black 
hole 1 and the other sources of energy (m 2 , q 2 )- 


E. Higher orders 

At higher orders, a proof of the extremum principle ([^ 
can be expected to proceed similarly to Sec. mmi so this 
subsection provides additional comments, rather than a 
full analysis. Through second order, as already noted in 
Sec. |IIID the mass m^ in (55) has the same form as the 
result (25) through first order: it is the same intrinsic 
function of the black hole’s area Ai and charge qi, inde¬ 
pendent of the other black holes (j ^ i). This result is 
perhaps unexpected, based on ([^, where m^ is a summa¬ 
tion that involves the other black holes (j ^ i). However, 
it could be anticipated physically, since the mass mi is 
specific to black hole i. Each mass mi contributes the 
the total energy. 


N 

E = mi 

i=l 


Eu 


(63) 


Our goal is to evaluate and extremize the energy function 
E{Ai,qi,rij). From the observations above, the sum of 


the masses in (63) can be anticipated to take a sim¬ 


ple form at higher orders. To evaluate E, it remains to 


consider the interaction energy, which we can view as the 
essential new task at each higher order, as a comparison 


of (26) and (57) illustrates; the steps will be similar to 


those following the evaluation of the area Ai in Sec. 


Through second order, our results of Sec. HID 


1 Sec. 

1—1 

1—1 

0 

HID 

show 


that that Aint = 0 at the energy extremum, if all charges 
have the same sign. This is seen by letting let qi = Emi 
for extremal charges, each with the same sign (±). This 
gives Qi = EM-i and 7) = ESi from ([4^ and (ImI), re¬ 


spectively. Then several terms cancel in the energy (57), 


which reduces to E = J2i shows that Ei„t = 0. 

It is promising to note that, at all orders, the follow¬ 
ing converse of the above result holds: if all charges are 
extremal and have the same sign, then the interaction 
energy vanishes, Ajnt = 0. This is most easily seen by 
substituting the extremal charge condition in the form 
0 into the interaction energy ([^. Since this is an ex¬ 
act result, it must hold at all orders, when all quantities 
are expanded in the inverse separations. The vanishing 
of Eint for extremal same-sign charges is an exact sup¬ 
porting result, and is a precise statement about energy, 
like the extremum principle. 


IV. ENERGY MINIMUM 

For sufficiently well separated bl ack ho les, as spe cified 
by the condition 17^ > 0 (see Secs. IIIC and HID), it is 


straightforward to show that our energy extremum, (29) 


and (60), is a minimum. To show this, we verify that the 


second derivatives of E are positive. Differentiating (27) 


and (58), holding constant the charges qi and separations 


gives 


d‘^E d'^m. 


dUi dmi 


dAi'^ dAi dAi 


(64) 


We then evaluate U^mi/dAi^ from (25) or (46), and we 


evaluate at the extremum {dmi/dAi = 0). This gives 
d^E 


Ui 


dA/^ 8v^Aj3/2 ’ 


(65) 


where Ui is to be evaluated with = |gi|. The condition 
Ui > t) ensures that (65) is positive, hence the energy 


extremum is a minimum. 

If each charge qi has the same sign, it is also straightfor¬ 
ward to verify that at the extremum, E = \Q\, where Q 
is the total charge. This is seen as follows. Let qi = Emi 
for extremal charges, each with the same sign (±). Then 


N 


N 


Q = = ±^m,. 


( 66 ) 


2=1 


2 = 1 


As already noted in Sec. |IIIE[ for extremal same-sign 
charges {qi = Emi), we also have 


N 


E = '^mi = EQ = \Q\. 


(67) 
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The energy minimum and its extremum value are both 
consistent with the BPS bound in the energy inequality 
[5] that is satisfied by the Hartle-Hawking static solution 
[2]. In contrast to [S], we have minimized the energy 
E{Ai,qi,rij) as a function of physical variables, without 
using supersymmetry. 

V. CONCLUSION 

Our main result in this paper is an extremum principle, 
which derives the extremal charge condition (jg^l = m^) 
for a set N black holes (all initially at rest and arbitrar¬ 
ily charged) by extremizing the total energy E{Ai, qi, Vij) 
with respect to the black hole horizon areas, at fixed 
charges and Euclidean separations. This principle is mo¬ 
tivated by the first law of black hole mechanics, and is 
valid if the black holes are sufficiently well separated, as 
specified by the condition Ui > 0, where Ui can be inter¬ 
preted as the effective gravitational potential experienced 
by black hole i. 

Our energy extremum is taken at fixed Euclidean sep¬ 
arations Tij, rather than fixed proper distances between 
horizons. This is natural, in the sense that the proper 
distances are well known to become infinite for extremal 
charges, while the separations remain finite. 

If all of the black hole charges have the same sign, then 
our extremum principle is a type of variational principle, 
and augments a list of existing principles [MS] that inter¬ 
pret static black holes (here, the static Hartle-Hawking 
solution) as extrema of total energy. Our results are also 
consistent with the supersymmetric BPS energy mini¬ 
mum {E = IQI) of [9|. Our derivation also shows how 
the corresponding substructure arises (rrii = |gi|) for each 
individual black hole. To our knowledge, our results pro¬ 
vide the first energy interpretation of the static Hartle- 
Hawking solution [2| . This interpretation requires neither 
the use of balanced forces, nor supersymmetry. 

It would be interesting to consider the regime of very 
small black hole separations, which is beyond the scope of 


this paper. This would probably require numerical meth¬ 
ods, and would involve an outermost common apparent 
horizon surrounding two nonextremal black holes that 
are sufficiently close to each other, similar to the case of 
uncharged black holes released at rest 0 [15] . The for¬ 
mation of a common apparent horizon has been studied 
numerically in the head-on collision of symmetric like- 
charged black holes m for charge-to-mass ratios in the 
range 0 < q/m < 0.98. 

In this paper, to locate a black hole’s apparent hori¬ 
zon analytically, we have done so perturbatively, which 
is a well-known feature of the Brill-Lindquist geometry 
[5] that we have employed. One might wonder if a dif¬ 
ferent geometry could be used instead; this appears to 
be unlikely. For example, in the conformally flat geom¬ 
etry found from Misner’s well-known method of images 
[M ITT] . each black hole’s apparent horizon is designed 
to be an exact coordinate sphere in the flat background 
space, which is analytically convenient. This is achieved 
by constructing the solution as an infinite series (sim¬ 
ilar to the method of images in electrostatics). How¬ 
ever, this approach does not permit black holes with ex¬ 
tremal charge, since the infinite series solution must sat¬ 
isfy a convergence condition. This condition is rather 
formal in general; for two symmetric black holes with 
opposite charges, it reduces to the statement that the 
black holes are nonextremal m- A similar restriction 
can be expected for same-sign charges. This suggests 
that the Brill-Lindquist geometry [5] is the unique fam¬ 
ily of nonextremal solutions that smoothly connects to 
the extremal Hartle-Hawking solution. 
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